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1 Introduction 



Tensions in several phenomenological models grew with experimental results on neu- 
trino/antineutrino oscillations at Short-Baseline (SBL) and with the recent, carefully 
recomputed, antineutrino fluxes from nuclear reactors. At a refurbished SBL CERN- 
PS facility an experiment aimed to address the open issues has been proposed p], 
based on the technology of imaging in ultra-pure cryogenic Liquid Argon (LAr). 

Motivated by this scenario a detailed study of the physics case was performed. 
We tackled specific physics models and we optimized the neutrino beam through a 
full simulation. Experimental aspects not fully covered by the LAr detection, i.e. the 
measurements of the lepton charge on event-by-event basis and their energy over a 
wide range, were also investigated. Indeed the muon leptons from Charged Current 
(CC) (anti-)neutrino interactions play an important role in disentangling different 
phenomenological scenarios provided their charge state is determined. Also, the study 
of muon appearance/disappearance can benefit of the large statistics of CC muon 
events from the primary neutrino beam. 

Results of our study are reported in detail in this proposal. We aim to design, 
construct and install two Spectrometers at "NEAR" and "FAR" sites of the SBL 
CERN-PS, compatible with the already proposed LAr detectors. Profiting of the 
large mass of the two Spectrometers their stand-alone performances have also been 
exploited. 

Some important practical constraints were assumed in order to draft the proposal 
on a conservative, manageable basis, and maintain it sustainable in terms of time- 
scale and cost. Well known technologies were considered as well as re-using parts of 
existing detectors (should they become available; if not it would imply an increase of 
the costs with no additional delay). 

The momentum and charge state measurements of muons in a wide energy range, 
from few hundreds MeV to several GeV, over a > 50 m 2 surface, is an extremely chal- 
lenging task if constrained by a 10 (and not 100) millions € budget for construction 
and installation. Running costs have to be kept at low level, too. 

The experiment is identified throughout the proposal with the acronym NESSiE 
(Neutrino Experiment with Spectrometers in Europe). 

In the next Section the relevant physics pleas for the Spectrometer proposal are 
summarized. In Section [3] an extensive discussion of the physics case and the possi- 
ble phenomenological models are described. In Section [4] details of a complete new 
simulation of the beam fluxes and horn designs are reported to work out a realistic 
situation and to be possibly used as contribution to the eventual design study group. 
In Section [5] the choice and design of the Spectrometers are discussed. After a brief 
reminder of the details of the Monte Carlo simulation for neutrino events (Section [6J 
the obtained physics performances are outlined in Section [7} Sections [8] and [9] deal 
with the technical definition of the mechanical structure and electrical setting-up for 



the magnets. Sections n0\ and 11 show the use of possible detectors either with 



a coarse resolution (~1 cm) or with a boosted one (~1 mm). Section 12 debates 
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about background levels to be taken into account for the data taking, described in 
Section 13 The last two sections report about CERN setting-up, schedule and costs. 
Finally conclusions are recapped. 



2 Aims and Landmarks 

The main physics topic of a new SLB CERN-PS experiment is to confirm or reject 
the so called neutrino anomalies and in case of a positive signal fully exploit the new 
physics connected with sterile neutrinos. The reported anomalies are the — > V e 
transitions observed by LSND j2], the v e disappearance claimed by the Gallium source 
experiments for the Gallium solar neutrino reactions [3] and the V e disappearance of 
reactor antineutrinos [I]. 

If such anomalies are not connected with new physics the LAr plus Spectrometer 
experiments will report null signal in the v e and V e spectra, both in the absolute 
flux and energy shape, ruling out the anomalies at exceeding a's. In fact as in some 
sterile neutrino models, discussed in the following, cancellations occur between v e (z/ e ) 
appearance, via — > V e transitions and v e (z/ e ) disappearance, a null result in the 
v e sector by LAr experiments would not be sufficient to fully exclude sterile neutrino 
models. In this case it will be crucial to precisely measure disappearance since it 
would be impossible for any model to fit both null results in the v e and the z/^ sectors. 

We note that the best limit available so far on disappearance in the Am 2 range 
under study, published by CDHS [5] , is based on the analysis of 3300 events at the Far 
detector, corresponding to just one day of data taking of the experiment we propose. 

We therefore clinch that in case of null signal the conjunct experiments will be 
able to rule out all the existing anomalies as connected to the presence of sterile 
neutrinos. 

If instead the anomalies are really hints of the existence of sterile neutrinos the 
experiments will be able to: 

• provide a direct signature of non-sterile neutrinos through the detection of Neu- 
tral Current (NC) disappearance. In fact while CC disappearance may occur 
via active neutrino oscillations^ NC disappearance can only occur when active 
neutrinos oscillate to sterile ones. It goes without saying that NC disappear- 
ance should be at the same rate of z/^ charged current disappearance rate as 
measured by a spectrometeiR 



a A //-neutrino oscillating to a r-neutrino will indeed disappear because at these energies r leptons 
cannot be created by neutrino interactions as they are below threshold. 

b It is amazing that the evidence of sterile neutrinos via the disappearance of NC neutrino inter- 
actions, might be made with the same neutrino beam used 40 year ago to detect for the first time 
the existence of NC neutrino interactions. 
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• Characterize the sterile neutrino model, namely to determine: 

— all the parameters of 3+1 models (3 parameters: U± e , U^, Am^) or 3+2 
models (7 parameters: U^e, U^, Amf 4 , {7 5e , U^, Am^, A s ); 

— the number of steriles from the number of Am 2 needed to fit the data that 
will correspond to — > v e , — > V e appearances and disappearance. 
Furthermore any manifestation of CP violation will be an evidence of at 
least 2 sterile neutrinos participating to oscillations. 

— CP violation in the sterile sector from the comparison of — > v e and 

— > v e transitions. 

Therefore in case sterile neutrinos indeed exist the experiments will fully exploit 
the very rich and exciting physics information related to them. 

From this concise discussion that will be more extensively addressed in the follow- 
ing, it clearly appears that a spectrometer capable of measuring the muon momenta 
would complement the LAr experiment by allowing to: 

• measure disappearance in the entire available momentum range. This is a 
key information in rejecting/observing the anomalies over the whole expected 
parameter space of sterile neutrino oscillations; 

• measure the neutrino flux at the Near detector, in the full muon momentum 
range, which is decisive to keep the systematic errors at the lowest possible 
values. 

The measurement of the sign of the muon charge would furthermore enable to: 

• separate from in the antineutrino beam where the contamination is 
relevant. Such a measure is crucial for a firm observation of any difference 
between — \ u e and — > V e transitions, a possible signature of CP violation; 

• and finally reduce the data taking period by a concurrent collection of both 
and events. 
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3 Physics 



3.1 Introduction 

The recent re-computation of antineutrino fluxes from nuclear reactors jl] leds to the 
so called "reactor neutrino anomaly" [HI [7], i.e. a ~ 3% deficit of the reactor antineu- 
trino rates measured at short baselines by past reactor experiments, reinforcing the 
case of sterile neutrinos. 

The latest fits of cosmological parameters are in favor of a number of neutrinos 
above 3 (see e.g. [H]), provided that masses of extra neutrinos are below ~1 eV/c 2 [§]• 

Besides reactor antineutrinos [7J [TU] and cosmology, sterile neutrinos were invoked 
to explain results from beam dump experiments, LSND [2], accelerator neutrinos 
(and antineutrinos) as measured by the MiniBooNE experiment [PI) [12] and source 
calibration data of the Gallium solar neutrino experiments [3] . Global fits to neutrino 
oscillation data [13] are performed by adding just one sterile neutrino to the three 
active, "3+1" models, or adding two sterile neutrinos, "3+2". The latter models (see 
Fig. [TJ provide a reasonable good fit to the data. In particular "3+2" models can 
introduce CP violation in the sterile sector, explaining the discrepancy between the 
V e appearance detected by the LSND and the lack of v e appearance as reported by 
MiniBooNE. "3+1" models cannot introduce CP violation and when compared to 
"3+2" they are disfavored at 97.2% C.L. according to [13J. A way-out to improve the 
goodness of fit of the "3+1" model is to introduce the CPT violation |14j . 

Global fits do not provide a clear unique solution emerging from the data. Fur- 
thermore it should be noted that tension still exists between appearance and disap- 
pearance data [13]; for this reason even "3+2" global fits are not fully satisfactory. 

In this Section we set the physics case of an accelerator neutrino experiment 
aiming at being the conclusive one on this topic. We recall the capabilities of a LAr 
detector and we focus on the physics opportunities that a spectrometer could add to 
a liquid argon target. We first discuss measures like the LSND/MiniBooNE evidence 
of — > V e transitions and the v ei V e disappearances, as well as measures like z/ M , 

disappearances. Then the possible signatures of sterile neutrinos, including NC 
disappearance, will be addressed and eventually the predictions based on some models 
will be exploited. 

It is worth noting that, should the anomalies be confirmed, a short-baseline accel- 
erator neutrino experiment would have an outstanding set of discoveries in its hands, 
namely: 

• the discovery of a new type of particles, the sterile neutrinos, interacting only 
via the gravitational force; 

• the proof that at least two sterile neutrinos exist; 

• the proof that the oscillations between active and sterile neutrinos violate CP. 



7 



3.2 — > v e oscillations 

The detection of — > V e oscillations has been claimed by the LSND [2] and the 
MiniBooNE [12J experiments. The two results may be in agreement with each other. 
The KARMEN experiment [15], which did not report any evidence of these transi- 
tions, limited the LSND signal parameter space. Furthermore MiniBooNE did not 
confirm the result in the neutrino sector [llj. 

A spectrometer would add very little to the detection of electrons corresponding 
to the signature of — > V e transitions, which can be very well measured by a Liquid 
Argon detector. Indeed the Icarus Collaboration already reported in their previous 
proposal pQ an excellent sensitivity to such transitions. 

Nevertheless a spectrometer can play a fundamental role in the determination of 

— > V e transitions by measuring the muon charges. That corresponds to a notable 
constraint: in the negative focussed beam, as discussed in Sect. [4], the rate of 
interactions is a sizable fraction of z/ M interactions. It turns out that it is mandatory 
to disentangle and rates, possibly on an event-by-event basis, at the Near 
detector, in order to reduce the systematic errors associated with the prediction of 
the Vfj, and fluxes in the Far detector. 

Considering that the Icarus detector would collect the LSND statistics in about 
10 days, it is evident that its sensitivity will be dominated by systematic errors and 
the information added by the Spectrometers looks very important, if not strictly 
mandatory to keep them as small as possible. 

3.3 u fl and V fl disappearance 

No experiment has so far reported evidence of or disappearance in the allowed 
Am 2 region for sterile neutrinos. The limits set by CDHS [5], MiniBooNE [16J and 
atmospheric neutrino experiments [T7] are among the most severe constraints on 
sterile neutrino oscillations. 

We stress here that the CDHS limit is based upon the analysis of 3300 neutrino 
events collected in the Far detector of a two-site detection setup at the CERN-PS 
neutrino beam. Such number of events corresponds to the statistics collected by our 
experiment in just one day's data taking. 

An improved measurement of (u^) disappearance could severely challenge the 
sterile global fits in case of null result or provide a spectacular confirmation in case 
of signal observation. 

The Up and z/ M transitions are the main physics topics that a Spectrometer exper- 
iment could address. The measurement of and z/ M spectra at the Near detector in 
the full momentum range is mandatory to constrain systematic errors (the and 
flux ratios at the Near and Far detectors are expected to be different, as discussed in 
Sect.g. 

In Sect. [7] the computed sensitivity of a spectrometer in measuring and z/ M 
disappearance is also discussed. It is important to note that with just 3 years' running 
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with the negative focussing beam the experiment could improve the existing limits 
on disappearance and provide a measurement of z/ M disappearance, so far never 
performed in the sterile Am 2 region. 



3.4 NC disappearance 

Sterile neutrinos come into play to accommodate the anomalies in a global phe- 
nomenological framework by introducing a fourth light neutrino. Since the number 
of active neutrinos is limited to three by the measurement of the Z width [TS], addi- 
tional light neutrinos cannot have electroweak couplings (sterile neutrinos). 

The detection of (z7 M ) disappearance would be extremely important since it 
would also open the door to a sensitive search for the disappearance of NC events, 
which is a direct signature of the existence of sterile neutrinos. Indeed none of the 
anomalies reported so far can be interpreted as a direct manifestation of sterile neu- 
trinos. 

For instance may oscillate to v T that cannot produce the associated r lepton 
being under threshold for the CC interaction at the energies of the beam. 

Instead NC interactions can "disappear" only if active neutrinos oscillate to sterile 
neutrinos. NC events, either v e or can be efficiently detected by the Liquid Argon 
detector. However the transition rate measured with NC events has to agree with 
the CC disappearance rate once the — > v T and — » v e contributions have 
been subtracted (these rates are anyway small at the L /E values of the present beam 
configuration). Indeed the NC disappearance is better measured by the double ratio: 

NC 

CC Near /i N 

JVC "J 

CCFar 

The double ratio is the most robust experimental quantity to detect NC disappear- 
ance, once CCN ear and CCp ar are precisely measured thanks to the Spectrometers, 
at the Near and Far locations, via the disentanglement of and contributions. 



3.5 Modelization 

Anomalies and computed sensitivities were (and are) usually addressed with em- 
pirical two neutrino oscillation formulas. Recently, the interest in sterile neutrinos 
has been greatly reinforced because, after the appearance of the reactor antineutrino 
anomaly the LSND/MiniBooNE anomaly and the Gallium source anomaly can be 
accommodated together with all other existing measurements of neutrino oscillations 
in a single model which incorporates two sterile neutrinos [13J. While this model 
can provide a good overall x 2 by fitting existing data, tension still exists between 
appearance and disappearance oscillation results. 

As introduced in Sect. 3.1| the main physics reason why the "3+2" models provide a 



better overall fit with respect to the "3+1" models is that they settle the experimental 
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conflict between the — > V e signal as detected by LSND and the lack of — > v e 
oscillations as reported by MiniBooNE thanks to the introduction of CP violation 
that requires at least two sterile neutrinos. 

In the following, the discovery potential of a LAr+NESSiE experiment will be 
discussed in the context of "3+2" models by selecting their best fit points in the 
parameter space. For completeness the "3+1" models will also be considered. 

3.5.1 "3+2" neutrino oscillations model 

In a short-baseline (SBL) accelerator experiment, where Am^ ~ ^ m 3i ~ 0> the 
relevant v e appearance probability is given by [15] : 



P u ^ Ve =4|C/ e4 | 2 |[/ At4 | 2 sin 2 04i + 4|[/ e5 | 2 |[/ M5 | 2 sin 2 51 

+ 8 \ U e JJ^U e <jU^\ sin 04i sin 05i cos(0 54 - 5) , (2) 

with the definitions 

(fa = , S = arg (u: 4 U^U e5 U; 5 ) . (3) 

where symbols have the usual meaning. Eq. ^ holds for neutrinos; for antineutrinos 
one has to replace 5 — > —5. 

The survival probability in the same SBL approximation is given by 




P Va ^ Va = 1 - 4 1 - \ \U ai \ 2 > j |L/ m | 2 sin 2 a -4\U a4 \ 2 \U a5 \ 2 sin 2 054 (4) 



where 0y is given in Eq. (j3]). 

The x 2 °f th e "3+2" model fit to the SBL oscillation data is shown in Fig. [I] It 
appears that there is more than one solution for the oscillation parameters. In the 
following, the four best fits indicated in [13] plus the best fit point of the preliminary 
updated analysis by Karagiorgi et al. [21] and the best fit by Giunti-Laveder [21] are 



considered (see Tab. 1 : ) 



3.6 Probabilities at the Far detector 

The oscillation probabilities P{y^ —> v e ), P(y e — > z/ e ) and Piy^ — > v^) are computed 
at a distance of 850 m from the proton target for the six best values reported in 

c It could be argued that the x 2 values reported in Tab. [I] are "too good". This happens because 
several data are basically not sensitive to sterile oscillations (i.e. Chooz data 23 ), and are well 
fitted under any hypothesis. 
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Figure 1: y 2 from global SBL data (upper panel) and from SBL reactor data alone (lower panel) 
for the 3+1 (blue) and 3+2 (red) scenarios. Dashed curves were computed using the old reactor 
antineutrino flux prediction [20], solid curves are for the new one J^]/. All undisplayed parameters 
are minimized over. The total number of data points is 137 (84) for the global (reactor) analysis. 
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0.201 
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0.1150 


0.101 


1.39 


114.4/130 


3) 


1.00 


0.133 


0.162 


1.60 


0.151 


0.078 


1.48 


114.4/130 


4) 


0.90 


0.123 


0.163 


6.30 


0.135 


0.091 


1.67 


115.0/130 


5) 


0.92 


0.14 


0.14 


26.60 


0.077 


0.15 


1.7 


182.6/192 


6) 


0.90 


0.158 


0.152 


1.61 


0.130 


0.078 


1.51 


91.6/100 



Table 1: Parameter values and x 2 a t the global best fit points for the four best "3+2" fit points of 
\13)j . table entries 1-4, the best fit of Karagiorgi et al. ]21j . entry (5); (Am 2 's in eV 2 ) and the best 
fit by Giunti-Laveder \241 , entry (6). 
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Tab. [T] (see Fig. [2]). Many interesting features occur within the interplay of appear- 
ance/disappearance either for or v e sectors, the most exquisite distinctiveness being 
that the scenario may be rather complicate. Therefore to disentangle that scenario 
we will need the best measurements of and v e as well as and V e in the widest 
available energy range. 

In any case we like to underline some attractive outcomes from these probability 
computations, primarily for the "3+2" models: 

• in the electron sector competition occurs between v e disappearance and v e ap- 
pearance through the — > v e transitions. This is the tricky way by which 
"3+2" may cancel out v e appearance in MiniBooNE. 

• The different results in v e and V e by MiniBooNE may be explained by the 
behavior of the — > V e transitions which take into account the CP violation 
introduced by the "3+2" models. 

• A low energy v e excess emerges. The peak is entirely due to the interference 
term of Eq. [2] and it would be a signature of the existence of not one, but two 
sterile neutrinos. 

It should be noted that the MiniBooNE low energy excess detected in the neu- 
trino run can be accommodated by the "3+2" models. However when disap- 
pearance data are also considered the peak of the interference term does not fit 
anymore the measure |13j . 

• A sizable disappearance probability as large as 15% at the oscillation maxi- 
mum is predicted below l-j-2 GeV, depending on the different best fits. 

3.7 Neutrino Rates at the Far detector 

Neutrino interaction rates are computed using the neutrino flux discussed in Sect. |4j 
the GENIE |25j cross sections and the above-mentioned probabilities. Event rates 
are normalized to 2 years' run of the positive focussing neutrino beam, with 30 kW 
proton beam power, and 3 years' run of negative focussing neutrino beam. A neutrino 
efficiency of 100% is considered, while energy resolution effects and systematic errors 
are not included in the plots^J 

By considering the six test points of Tab. [T] event rate spectra are displayed in 
Fig. [3] for v e . Results without oscillations, with v e disappearance only and with both 
v e disappearance and appearance, generated by — >■ v e transitions, are shown. 

The number of expected events are displayed in Fig. Ill while event rates in the 
energy range 0.2 GeV < E v < 2 GeV are reported in Tab.[2]and Tab. [3j for and 
u e , respectively. The integral effect is not negligible (up to 6%) and the distinctive 
spectral signature is clearly detectable. 

d Results with full simulation, including energy resolution effects and systematic errors, are re- 
ported in Sect. [7] 
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Figure 2: Oscillation probabilities computed at a baseline of 850 m for the four best fit points of 
"3+2" fl3j . the best fit of l21\/ and the best fit of pffl . P(u^ -t v e ) on the left (neutrinos in blue, 
antineutrinos in magenta), P(y e — > v e ) at center and Piy^ —¥ v^) on the right. 
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Figure 3: Electron neutrino spectra for the six test points described in the text. Blue lines cor- 
responds to no oscillations, magenta to disappearance only, black points with statistical error bars 
corresponds to disappearance plus appearance events. 
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Figure 4: i/„ event rates (left) where blue lines correspond to no-oscillation and black points with 
statistical error bars (not visible) include disappearance. Difference between rates estimated with and 
without oscillation (right). The six rows refer to the test points discussed in the text. 
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Fit point No-Osc. Disappearance 

1) 605792 569060 

2) 605792 576102 

3) 605792 560108 

4) 605792 564136 

5) 605792 554629 

6) 605792 567584 

Table 2: events computed for 0.2 < E v < 2 GcV and two years' run at 30 kW, for the 6 test 
points described in the text. 



Fit point No-Osc. Disapp. Disapp. + App. 



1) 1424 1349 2065 

2) 1424 1342 1512 

3) 1424 1306 1788 

4) 1424 1326 2018 

5) 1424 1349 1977 

6) 1424 1283 1835 



Table 3: v e events with 0.3 < E v < 2 GeV computed for two years' run at 30 kW, for the 6 test 
points described in the text. 
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3.8 Antineutrino Rates at the Far detector 

If CPT holds, data taken with antineutrino beams should provide identical v e and 
disappearance rates. If CP is violated, as predicted by all the "3+2" best fit points, 
V e appearance rate results to be rather high and the signal could not be missed. As 
an example we display in Table [4] event rates as predicted by the "3+2" best fit and 
in Figure [5] the corresponding signal plots. 

Section [7] quantitatively discusses the possibility of measuring both and z/ M 
disappearance rates with the NESSiE Spectrometer in the antineutrino run. 



events 


No-Osc. 


Disappearance 


Disapp. + App. 


v e (0.3-2 GeV) 


1134 


1078 


1475 


(0.3-2 GeV) 


312439 


290791 





Table 4: Event numbers in antineutrino mode, best fit (1), 3 years' of data taking 




E (GeV) E(GeV) E (GeV] 



Figure 5: Left panel: v e +V e events in the antineutrino beam, computed for the first test point of 
Table 1, blue line corresponds to non- oscillated events, magenta line to disappearance only and points 
with errors correspond to oscillated events. Central panel: +V^ L events, blue line corresponds to 
non- os ciliated events, points with errors (not visible) to oscillated events. Right panel: Ua +V ^ events 
are shown as expected minus measured. 



3.8.1 "3+1" neutrino oscillations model 

In the so called "3+1" model the flavor neutrino basis includes three active neutrinos 
u e , Up, v T and a sterile neutrino v s . The effective flavor transition and survival prob- 
abilities in short-baseline (SBL) experiments are given using the standard notation 
by 

i? s f V, = sin 2 sin 2 (^^) (a*0), (5) 
^fV) = l-sin 2 2^sin 2 f^|^) , (6) 
for a, (3 = e, /i, r, s, with 
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siYi 2 2$ a p = A\U a4 \ 2 \U pA \\ (7) 
sm 2 2^ Q = 4|f/ Q4 | 2 (l-|£44| 2 ) • (8) 

The key features of this model are: 

1. All effective SBL oscillation probabilities depend only on the absolute value of 
the largest squared-mass difference Am^. 

2. All oscillation channels are open, each one with its own oscillation amplitude. 

3. The oscillation amplitudes depend only on the absolute values of the elements 
in the fourth column of the mixing matrix, i.e. on three real numbers with sum 
less than unity, since the unitarity of the mixing matrix implies Yl a \ U a 4\ 2 = 1 

4. CP violation cannot be observed in SBL oscillation experiments, even if the 
mixing matrix contains CP-violation phases. In other words, neutrinos and 
antineutrinos have the same effective SBL oscillation probabilities. 

The global fit of all data in "3+1" scheme yields the best-fit values of the oscillation 
parameters listed in Tab. |5j 



3+1 


V 2 . 
A min 

NDF 


100.2 
104 


GoF 


59% 


Am^ [eV 2 ] 
|t/ e4 | 2 


0.89 


0.025 




0.023 


NDFp G 


24.1 
2 


PGoF 


6 x 10~ 6 



Table 5: Values ofx 2 , number of degrees of freedom (NDF), goodness- of -fit (GoF) and best-fit 
values of the mixing parameters obtained in 3+1 fits of short-baseline oscillation data. The last 
three lines give the results of the parameter goodness- of- fit test, A\pq, number of degrees of freedom 
(NDFpq) and parameter goodness- of- fit (PGoF). 

Figures [6] and [7] show the allowed regions in the sin 2 2 i # e/J -Am| 1 , sin 2 2$ ee -Am 2 1 
and sin 2 planes, respectively, together with the marginal A% 2 's for Am 2 1; 

sin 2 2$ eM , sin 2 2 r d ee and sin 2 2i9 M/Lt . 

The proposed NESSiE experiment aims at exploring these regions. 
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10" 4 10" 3 10" 2 10" 1 2 4 6 8 10 
sin 2 2$ e[i Ax 2 



Figure 6: Allowed regions in the sin 2 2z9 eM -Am^ plane and marginal A% 2 's for sin 2 2z9 eM arte? 
Am|j obtained from the global fit of all the considered data in 3+1 schemes. The best-fit point 
corresponding to Xmin * s indicated by a cross. The isolated dark-blue dash-dotted contours enclose 
the regions allowed at 3<7 by the analysis of appearance data (the — > D e data of the LSND, 
KARMEN and MiniBooNE experiments and the -> v e data of the NOMAD JMj and MiniBooNE 
experiments). 
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Figure 7: Allowed regions in the sin 2 2t? ee -Am| 1 and sin 2 2i? AIM -Am 2 1 planes and marginal A\ 2 's 
for sin 2 2i9 ee and sin 2 2$^ obtained from the global fit of all the considered data in 3+1 schemes. 
The best-fit point corresponding to Xmin * s indicated by a cross. The line types and color have the 
same meaning as in Fig. [6| The isolated dark-blue dash-dotted lines are the 3er exclusion curves 
obtained from reactor neutrino data in the left plot and from CDHS and atmospheric neutrino data 
in the right plot. The isolated dark-red long-dashed lines delimit the region allowed at 99% C.L. by 
the Gallium anomaly. 
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3.9 3+1 and CPT violation model 



The only implementation among 3+1 models able to fit global data is the 3+1 and 
CPT violation model of Giunti-Laveder [33 14^| The model was inspired by the 
analysis of the electron neutrino data of the Gallium radioactive source experiments 
and the electron antineutrino data of the Bugey |22j and Chooz [23] reactor experi- 
ments in terms of neutrino oscillations allowing for a CPT-violating difference of the 
squared-masses and mixings of neutrinos and antineutrinos. 

It was found that the discrepancy between the disappearance of electron neutrinos 
indicated by the data of the Gallium radioactive source experiments and the limits on 
the disappearance of electron antineutrinos given by the data of reactor experiments 
reveal a positive CPT-violating asymmetry of the effective neutrino and antineutrino 
mixing angles. If there is a violation of the CPT symmetry, it is possible that the 
effective parameters governing neutrino and antineutrino oscillations are different. 
From a phenomenological point of view, it is interesting to consider the neutrino 
and antineutrino sectors independently, especially in view of the experimental tests 
considered in this proposal. 

The parameters of the model are reported in Tab. [6j 













-r '2 

Am 41 










1.92 0.275 


0.0 


0.47 0.061 


\ 0.886 



Table 6: Best fit parameters of the 3+1 and CPT violation model 

In this scenario, neutrinos undergo u e — > u e transitions only, see Fig. [8j while 
antineutrinos have a much richer phenomenology [28| [29] . 

The v e spectra computed for two years' run at 30 kW are reported in Fig. [9j 
According to 3+1 and CPT violation model v e disappearance should be clearly de- 
tectable 1127 events detected against a prediction of 1424 for 0.3 < E u < 2 GeV 
(Tab.0). 



events 




No-Osc. 


Disapp. 


Disapp. + App. 


v mode 


p e (0.3-2 GeV) 


1424 


1127 


1127 


V mode 


v e (0.3-2 GeV) 


1260 


1254 


1699 


V mode 


(0.3-2 GeV) 


329328 


277335 





Table 7: Events expectation for the 3+1 and CPT violation model 



c Indeed also"Non Standard Neutrino Interactions" or quantum decoherence have been pro- 
posed [27]. 
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Figure 8: Oscillation probabilities computed for the 3+1 and CPT violation model. P{v^ — > v^) 
and P{v^ — > v e ) are not displayed because they are predicted to be null by the model. 




Ev (GeV) Ev (GeV) 

Figure 9: Top panel: v e event rates in the negative focussing beam (left) and V e event rates in the 
positive focussing beam (right) computed assuming the 3+1 and CPT violation model. N on- oscillated 
rates are displayed as a blue histogram. Bottom panel: event rates in the positive focussing beam: 
on the right absolute event rates computed with (point with errors) and without (blue histogram) 
oscillations following the 3+1 and CPT violation model; on the left the difference between expected 
and measured events is shown. 
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4 The Neutrino Beam 



The NESSiE detector is planned for exposure to the CERN-PS neutrino beam-line, 
originally used by the BEBC/PS180 collaboration [30] and re-considered later by 
the I216/P311 [31J proposal. The baseline setup used by BEBC/PS180 consisted 
in a 80 cm long, 6 mm diameter beryllium-oxide target followed by a single pulsed 
magnetic horn operated at 120 kA. The PS can deliver 3 • 10 13 protons per cycle 
at 19.2 GeV kinetic energy in the form of 8 bunches of about 60 ns in a window of 
2.1 fis, integrating about 1.25 • 10 20 protons on target (p.o.t.)/year under reasonable 
assumption^} The existing decay tunnel, which has a cross section of 3.5 x 2.8 m 2 
for the first 25 m of length and 5.0 x 2.8 m 2 for the remaining 20 m, is followed by a 4 
m thick iron shield and 65 m of earth. With respect to the original configuration, the 
target and the horn must be redesigned and reconstructed due to the present level 
of radioactivity while the proton beam line magnets and supplies could be recovered 



(see Sect. 4.3) 



4.1 Beam simulation 

The BEBC/PS180 fluxes were reproduced by I216/P311 in their Letter of Intent [32] 
using a simulation based on GEANT3 and GFLUKA. The comparison was updated 
using a simulation adapted from the one used for the CNGS beam for their pro- 
posal [31]. For this memorandum a GEANT4 and FLUKA [33] based simulation 
has been developed to profit of a modern programming framework and to investigate 
possible improvements in the beam performance. 

The generation of proton-target interactions is done with FLUKA-2008 while 
GEANT4 is used for tracking in the magnetic field and the materials and for the 
treatment of meson decays. The simulation program, thoroughly described in [34] . 
has been modified to take finite-distance effects into account, which are particularly 
important due to the short baselines involved (127 and 885 m measured from the 
target to the beginning of the LAr detector). Neutrinos crossing the LAr and Spec- 
trometer volumes are directly scored using a full simulation avoiding any weighting 
approach)^} A sample of 10 7 simulated p.o.t. was produced and has been used in the 
following. 

In order to benchmark the GEANT4 simulation the existing setup used for the 
BEBC/PS180 experiment has been reproduced [35]. The layout of the considered 
volumes is shown in Fig. [To] 

The spectrum shape of the z/^ is in good agreement with that calculated by the 
I216/P311 experiment whereas the obtained normalization is instead 27% lower (Fig. 



11). Investigations are ongoing to understand this difference, in particular the im- 



plementation of the geometry and the hadro-production models are under checking. 



f 180 days' run per year allocating one third of the protons are assumed, with present perfor- 
mances. 

s This technique is relatively CPU-consuming but nevertheless affordable. 
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Figure 10: Layout of the focussing system and Near detector station in the GEANT4 simulation. 



The flux reduction resulting from negative-focussing, which amounts to about 40%, 
is well reproduced. 

Neutrino fluxes at the Far and Near Spectrometers in positive and negative- 



focussing are shown in Fig. 12 Accounting for the Spectrometer geometries and 
locations we obtain a ratio u/p.o.t. of about 10~ 2 in the Near station and a further 
reduction of a factor of about 20 in the Far location. The fractions of (z/^, u e , u e ) 
in the Near Spectrometer are (92.3, 6.6, 0.87, 0.26)% in positive-focussing mode and 
(12.2, 86.6, 0.49, 0.75)% in negative-focussing mode. In the Far detector numbers are 
very close to the previous ones. The electron neutrino contamination is in agreeement 
with the estimates of I216/P311. 

The energy dependence of the electron and wrong-sign muon contamination is 



shown in Fig. 13 It can be noticed that the v e contamination in the energy region 
between 500 MeV and 1 GeV is very low, approaching a minimum of 0.4%, to the 
benefit of the v e appearance search. The contamination of in the beam is 
quite significant especially at high energy where the spectrometer charge separation 
becomes thus very important. 

In particular, after folding the spectrum with the CC cross section (Fig. 14), one 
can see that the negative-focussing beam at high energy yields a comparable mixture 
of fi~ and fi + . The Spectrometer charge separation allows studying the behavior of 
both CP states at the same time. 

The distribution of the impact points of the in the Near Spectrometer is shown 
in the upper left plot of Fig. 15 The center of the beam is displaced in the bottom 
direction of 1 m and the radial beam profile can be well fitted with a Gaussian having 
a cr of about 4 m (Fig. 15 upper right and lower left plots) 

The spectrum of at the Near detector peaks at lower energies due to the fact 
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Figure 11: Comparison of the present simulation (dots) with those of BEBC (dashed line) and 
I216/P311 (continuous line) for positive-focussing mode (left) and negative- focus sing mode (right). 
G4 stands for the new simulation obtained with GEANT4 which introduces a 1.27 reduction factor 
with respect to the old simulation. Fluxes are referred a distance of 825 m from the target (BEBC 
site). The target starts at 35 cm upstream of the BEBC/PS180 horn, which is operated at 120 kA. 
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Figure 12: fluxes (black) and fluxes (red) at the Near (left column) and Far (right column) 
Spectrometers in positive-focussing (upper row) and negative-focussing (lower row) modes. 
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Figure 13: Left: v e jv^ ratio in positive-focussing mode and v^jv^ ratio in negative-focussing mode 
in the Near Spectrometer as a function of the neutrino energy. Right: Cv/zAj ratio in positive- 
focussing mode and v^jv^ ratio in negative-focussing mode in the Near Spectrometer as a function 
of the neutrino energy. 



Event rates. M=0.6 ton, LAr, 2.5e20 pot, L=0.85 Km 




Event rates. M=0.6 ton, LAr, 2.5e20 pot, L=0.85 Km 




Figure 14: Comparison of v^f c and v^ c event rates for the positive (left) and negative-focussing 
(right) beams. 
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Figure 15: Up left. Distribution of the neutrino impact points in the transverse plane, X vs Y . 
Up right: X projection. Bottom left: Y projection. The black histogram represents the v^, the red 
histogram the v^. Bottom right: E v vs V ' X 2 + Y 2 . Results refer to positive-focussing mode. 



that many interacting neutrinos are off-axis and then tend to have a lower mean 
energy (Fig 16, left). Nevertheless, restricting to the central region, i.e. the region 
that subtends an angle similar to that of the Far detector, the shapes tend to get 
closer. This selection allows smaller corrections in the Near/Far ratio thus decreasing 
the systematics. The correlation between the energy and the radius of the neutrino 
impact point is shown the lower right plot of Fig. 15 The shape of the u e spectrum 
on the other hand is more similar in the Near and Far locations (Fig 16, right) due 
to the predominantly 3-body decay origin. 



4.2 New horn designs 

By keeping the basic geometry of the beam (size of the target and horn hall, decay 
tunnel, shielding) unchanged, studies of several modified horns and targets in addition 
to the BEBC original setup are undergoing, an optimized design of new and possibly 
improved beam optics being pursued. Our aim is to move towards a high intensity 

flux in order to reduce the statistical errors, peaking at {E„ ) ~ 1 GeV to match 
the Am 2 ~ 1 eV 2 region. For the reduction of systematics in the — > v e appearance 
channel the high energy tail above ~ 2.5 GeV should be reduced to suppress tt° 
production and the intrinsic v e contamination should be kept small. 

Preliminary studies performed with bi-parabolic horns d la NuMI [36J indicate 
room for improvement. In particular in the last decades the feasibility of horns pulsed 
at high currents of order 300 kA has been demonstrated. The fluxes obtained with a 
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Ev (GeV) E, (GeV) 

Figure 16: Comparison of the v„ (left) and v e (right) fluxes in the Near and Far LAr detector in 
positive-focussing mode. G4 stands for the new simulation obtained with GEANT4 which introduces 
a 1.27 reduction factor with respect to the old simulation. 



new bi-parabolic horn pulsed at 246 kA, are shown in Fig. 17 An overall flux gain is 
obtained, particularly at low energy. This would be particularly useful to extend the 
sensitivity towards smaller Am 2 . 

Since the interval of the interesting Am 2 regions is large it would be desirable to get 
a configuration capable of scanning different neutrino energy regions thus adapting 
to the physics indications coming from the data taking. Despite the fact that the 
overall features of the neutrino beam are dictated by the proton energy, an effective 
way for varying the mean energy consists in using a tunable target-horn distance. 

While keeping the shape of the original horn unchanged, a scan in the scatter plane 
(current (i) vs target longitudinal position (z tar g)) was performed to study the effects 
of the fluxes in terms of normalization, energy distribution and contaminations. 
Results are shown in Fig. 18 for positive-focussing at the Far location. Plots show how 
the integral flux is increased going to higher currents; pushing the target upstream 
(downstream) corresponds to decrease (increase) the contamination and to probe 
the high(low)-energy region. Two configurations yielding very different spectra are 



shown in Fig. 19 It is interesting to notice that these fluxes are obtained using the 



same current of 300 kA just varying the position for the target. 



4.3 Beam reactivation 

Preliminary evaluations for a renovated TT7 PS neutrino beam line have been al- 
ready completed at CERN [37J. The TT7 transfer line, the target chamber and the 
decay tunnel are in good shape and available for the installation of the proton beam 
line, target and magnetic horn. The main dipoles, quadrupoles, correction dipoles 
and possibly the transformer for the magnetic horn can be recuperated, reducing 
significantly the cost and time schedule. The target and secondary beam focussing 
design can profit of the CNGS experience as well as monitoring systems, primary 
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Figure 17: New horn fluxes compared to BEBC and I216/P311. The meaning of the symbols is as 
in Fig. in particular G4 stands for the new simulation obtained with GEANT4 which introduces 
a 1.27 reduction factor with respect to the old simulation. 

beam steering and target alignment |38j . 

The study dating back to 1999 [31] estimated the time required to be approxi- 
mately 2 years for a total cost of about 4.2 MSF detailed as follows: 

• power converters of TT1 and TT7 magnets: consolidation and lower level elec- 
tronics (1.1 MSF); 

• civil engineering, mostly new housing for converters (0.5 MSF); 

• removal of 400 m 3 radioactive waste material in TT7, provisions and installa- 
tions for radiation protections, access control (< 0.5 MSF); 

• beam line installation, vacuum chamber, general mechanics (0.4 MSF); 

• beam monitoring instrumentation (0.4 MSF); 

• new target and horn (0.4 MSF); 

• new pillars and platform in the Near experimental hall (50 KSF). 
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Figure 18: Dependence of the integral flux (top), the mean energy (middle) and the 
contamination (bottom) from (I, z tar g) with a test bi-parabolic horn. Results are shown for positive- 
focussing at the Far location. 
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Figure 19: Two extreme low (red) and high-energy (blue) configuration fluxes arising from the scan 
in the (I, z tar g) plane. The current is 300 kA and the position of the beginning of the target with 
respect to the beginning of the horn is +32 and -58 cm respectively. The fluxes refer to a distance 
of 825 m. 
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5 Spectrometer Design Studies 



The main purpose of a spectrometer placed downstream of the target section is to 
provide charge and momentum reconstruction of muons escaping from LAr detection. 
This choice would provide a double benefit with respect to a LAr detector running 
in stand-alone mode. Firstly a precise muon momentum reconstruction would allow 
a good kinematical closure of CC events occurred in LAr in particular in the high 
energy tail where the muon momentum resolution is poorer. Secondly muon charge 
separation would allow us to disentangle and P M disappearance channels in partic- 
ular in the negative-focussing option where wrong-sign contamination is larger. This 
would allow tackling CP and CPT violating scenarios in an unambiguous way. 

In addition, a proper mass-granularity combination would allow a coarse recon- 
struction of vz. events occurring within the Spectrometer itself and provide a way 
to use the Spectrometer in stand-alone mode, too. 

We characterize the detector by evaluating the physics performances in terms of 
CC disappearance and the sensitiveness to (Am 2 , sin 2 29) values predicted by the 
models described in Sect. [3j We recall that with a suitable detector the old CDHS 
disappearance limit will be tested in just one day of data taking. 

Moreover we assumed some "a priori" constraints by defining a realistic, conser- 
vative, relatively inexpensive apparatus. 



5.1 Magnetic field in Iron 

A magnetized iron spectrometer was chosen as baseline option. Its design was such to 
match the required sensitivity of the experiment with respect to the sterile neutrino 
models of Sect. [3} The basic parameters to be tuned are the transverse widtb|^J the 
longitudinal dimensions, the iron slab thickness and the tracking detector resolution. 
The transverse and longitudinal dimensions constrain the detector acceptance and 
hence the detectable number of events which are directly related to the and 
1 mixing parameters (see in particular Table [l] of Sect. [3). On the other hand the 
target segmentation sets bounds on muon momentum and charge reconstruction and 
are directly related to Am| x and Am 41 . 

The transverse width has to be large enough in order to maximize the detection 
of muons escaping from LAr. Monte Carlo simulation has shown that a transverse 
size matching the LAr acceptance (~ 8 x 5 m 2 in the Far site (FD) provides a good 
performance while still keeping the detector size at a feasible level (doubling the 
dimensions would improve the sensitivity to by 5%). 



Fig. [20] shows the muon impact point distribution along the x axis at FD. 
Due to the relatively low energy spectrum of the PS neutrino beam the Spec- 
trometer longitudinal size has to be large enough in order to allow energy measure- 



h Here and in the following the transverse directions are defined by the horizontal x and the 
vertical y axis while the z axis runs in the beam direction. 
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Figure 20: Distribution of muon impact points along the x axis, at the Far site (FD). The dip at 
the center of the distribution and the sharp tails on the sides are due to the convolution with the LAr 
detector acceptance. The shaded area corresponds to the 78% of the total one. 
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ment by range in a wide energy interval (see later on in this Section). Simulations 
have shown that 2 m of Iron contain ~90% of muons in positive-focussing (~85% in 
negative-focussing). 5 cm thick iron slabs would allow a spectrum bin size of ~60-=-100 
MeV/ a/12 ~ 20 -j- 30 MeV for perpendicularly impinging muons. Momentum recon- 
struction for passing-through muons (E^ > 3 GeV/c) can be performed exploiting 
the track bending in the magnetic field. A detector resolution of ~1 cm provides a 
(Jp/p resolution ranging from 20% at 3 GeV/c up to 30% at 10 GeV/c (see Fig. 21 ). 




10" 1 2x10" 1 3x10" 1 1 2 3456789 10 20 

Muon Momentum (GeV) 

Figure 21: Muon momentum resolution as calculated for several Spectrometer configurations, with 



a magnetic field B = 1.5 T . 



Muon charge assignment is performed measuring the direction of track curvature. 
The experimental sensitivity to Am 2 scales linearly with the neutrino energy, E u , and 
only with \/N events . Therefore the possibility to explore CPT-violating disappear- 
ance at Am 2 values below ~1 eV 2 in a baseline of ~ 850 m relies on the capability 
to identify the muon charge with good accuracy down to E u of about few hundreds 
of MeV and even lower if we consider the muon residual energy reaching the Spec- 
trometer (see Fig. 22). However in this momentum region one has to cope with the 
severe limits imposed by Multiple Scattering in iron which completely dominates the 
charge identification capability. 

Since for stopping muons the path-length L in the B ■ L product is fixed by range 
B has to be as large as possible (B ~ 1.5 T in iron). Once the magnetic field 
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Figure 22: Energy spectrum of fi + intercepted by the Spectrometer in the FD in negative-focussing 
for 3.75 p.o.t. The shaded region is the subsample of muons from events such that 1.27 x 1 eV 2 x 
L/Ec > it/2. 



35 



is fixed the momentum resolution u p /p and the charge mis-identification rj (defined 
as the fraction of muon tracks whose charge is wrongly assigned) are limited at low 
energies by Multiple Coulomb Scattering (MCS). In order to understand qualitatively 
this behavior let us define 77 as the fraction of muon tracks whose charge is wrongly 
assigned. 77 is related - in the Gaussian approximation of the Moliere distribution- to 
the momentum resolution by 



77 = -erfc 
' 2 



(9) 



where erfc is the complementary error function. The capability to separate charge- 
conjugated oscillation patterns depends on the available statistics and on the value 
of 77: the lower the number of events the lower 77 has to be in order to assess a sepa- 
ration at a given significance level. For instance neglecting systematic and statistical 
uncertainties on 77 a vjv separation at 10% level would require 77 ~ 30% in energy 
bins with 10000 events and 77 ~ 5% with 1000 events. 

The momentum resolution (or equivalently rj) is the combination of two terms, 
the first one related to MCS and the other to measurement errors: 




V V W J Mas V P 



(10) 



The two terms can be calculated e.g. for the case of uniform detector spacing as 
in [53 H01 HE]- The first term is almost independent of the number of measurement 
points along the trajectory and it is expressed as 

c„\ / 1.5T\ / lm^ ' " 



where L is the muon path-length in iron and is actually limited by the muon range 
and the spectrometer size at low and high momenta respectively. This term sets the 
lower irreducible limit that can be obtained by a measurement of this kind (note that 
77 ~ 5% roughly corresponds to cr p /p — 50%). 

The second term depends on the detector resolution and on the number of mea- 
surements (or equivalently on the iron slab thickness, A) and, in principle, can be 
decreased by changing the detector sampling and/or space resolution. In this case, 
for L » 5A, 



<\ * 13% (^a (ay 2 r^r r±n 



lGeV/cy V5cmy \ L J \ B J VlcnJ 



(12) 



with the same considerations on L as for equation 11 



Curves reported in Fig. 21 and Fig. 23 show the momentum resolution and the 77 



dependence on the muon momentum for various choices of A and adet (for orthogonally 
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Figure 23: Muon charge mis-identification rj as calculated for several Spectrometer configurations, 
with a magnetic field B = 1.5 T . 
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impinging muons). It is apparent that at best a mis-identification rj below 5% can be 
obtained only for muons with momenta above ~ 500 MeV/c. 

An option to improve charge identification capability below 500 MeV is to equip 
the Spectrometer with a magnetic field in air just in front of the first iron slab. This 
possibility is discussed in the next Section. 

From now on we will assume as baseline option for the FD site a dipolar magnetic 
Spectrometer made of two arms separated by 1 m air gap. Each arm consists of 21 
iron slabs 585x875 cm 2 wide and 5 cm thick interleaved with 2 cm gaps hosting 20 
Resistive Plate Chambers (RPC) detectors with ~1 cm tracking resolution. The Near 
detector (ND) Spectrometer is a downsized replica of the FD one, the only difference 
being in fact the transverse dimensions (351.5x625 cm 2 ). 

A calculation of the magnetic field map was performed with the COMSOL code 



Fig. 24 shows the distribution of the magnetic field in the 21+21 iron layers (for con- 
venience we show in the figure the entire magnetic system, i.e. in iron and in air, the 
latter to be discussed later on in the Section). The fringe field is 70 -j- 100 G at 10 cm 
distance from the iron slab edge (Fig. [25] ) . 

The simulation was performed assuming two symmetric coils wrapped around the 
top and bottom flux return path (following the same concept applied in the OPERA 
magnet [45J). The turns (36) are made of aluminum bars 27 x 19 mm 2 with inner 
water cooling. The current density is 3 A /mm 2 for a total current of 1255 A. The 
total resistance is 39 mOhm per coil, the voltage is 49 V and the power 61 kW. The 
conductor cross section can be increased to reduce the power dissipation. 

5.2 Magnetic field in Air 

For low momentum muons the effect of Multiple Scattering in iron is comparable 
to the magnetic bending and therefore the charge mis-identification increases (see 



Fig. 23 ) For muon momenta < 1 GeV/ c the charge measurement can be performed 



by means of a magnetic field in air. In Fig. 26 the displacement expected in the 
bending plane is shown for muons crossing a magnetized air volume of 30 cm depth. 
A uniform magnetic field oriented along the y axis (the bending plane is the z, x one) 
is assumed. In the left plot the shift is shown as a function of the muon momentum 
for some values of the magnetic field in the 0.1 — 0.4 T range. In the right panel 
the spatial displacement in the bending plane estimated for muons of 0.5 GeV in a 
magnetic field of 0.3 T as a function of the incoming angles is plotted. 

A simulation of the magnetic field in air was realized based on a coil wound on a 
large conductor (54 mm x 19 mm) with 170 turns distributed along the spectrometer 
height. A uniform magnetic field of 0.25 T along the y axis is obtained with a coil 



current density of 8 A/mm 2 (Fig. 24, 7000 A current). The fringe field at 10 cm 
distance from the edge (Fig 27) is 70 -j- 200 G. The total resistance is 0.1 Q, the 
voltage is 700 V and the power 4.9 MW. This set of parameters shall be intended 
as preliminary and needs to be optimized for reducing the electrical power and for 
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Figure 24: Simulated magnetic field distribution in the Iron Spectrometer (yellow 



regions) and in the air volume (light blue region on the left side, see Sect. 5.2). The 
horizontal axis is in meters; the color bar is in Tesla. The iron field is generated by 
two coils wrapped around the top and bottom iron flux return (only the upper coil is 
shown in the figure). 
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Figure 25: Magnetic field distribution in the 21+21 iron slabs and in the air volume 
upstream of the first iron slab. The vertical axis indicates the magnetic field in Tesla. 




muon momentum (GeV/c) angle in x,y plane O 

Figure 26: Left: Spatial displacement of muon trajectory in the bending plane as a function of their 
momentum after crossing an air magnetic field volume of 30 cm thickness. Different field values in 
the 0.1-0. 4 T range are considered. The incoming muon direction is perpendicular to the detector 
planes. Right: the displacement of 0.5 GeV muon tracks in the x,z view (B = 0.3 T) as a function 
of the initial muon angle (the color bar scale is in cm). 
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splitting the single large coil into a number of smaller (in height) and more manageable 
coils. 




Figure 27: Simulated magnetic field distribution in the Iron Spectrometer (yellow 
regions) and in the air volume (light blue region on the left side). The horizontal axis 
is in meters. The color bar scale is in Tesla. The vertical pairs of coils used to generate 
the field in air are shown (on the left of the picture). 



As discussed in Sec. 11 charge measurement for low momentum muons can be 
performed by designing a precision tracker detector of ~ 1 mm spatial resolution 
installed in the z, x plane. The reconstruction in the z, y view would be not required 
since the momentum could be measured by the muon range in the Iron Spectrometer. 
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6 Monte Carlo Detector Simulation and Recon- 
struction 



The present proposal has been extensively developed using full-detail programming 
and up-to-date packages to obtain precise understanding of acceptances, resolutions 
and physics output. Although not all possible options have been studied, a rather ex- 
haustive list of different magnet configurations and detector designs has been adopted 
as benchmark for further studies once the detector structure will be finalized. 

6.1 Simulation 

The aim of the simulation of the apparatus is to help the design studies reported here 
to understand the main features of the proposed experiment. The simulated detector 
consists of a ND and FD part, both made of an Liquid Argon target followed by a 
magnetic Spectrometer. The relative position and the dimensions of the LAr target 
have been kept fixed with respect to the Spectrometer for all the studies described 
below. 

Muon neutrino and antineutrino fluxes for positive and negative beam polarity 
were assumed as those reported in the LAr proposal [1] since the beam analysis was 
still in progress. In the simulation the beam has 1 mrad tilt with respect to the 
horizontal. For the time being no angular dependence of the neutrino fluxes has been 
considered. For the future we plan to take profit of the full simulation of the secondary 
beam line reported in the previous Sect. [4] in order to improve the simulation of the 
neutrino beam at both the Near and Far detector positions. 

Neutrino interactions are generated in the Argon target using GENIE [25] with 
standard parameters and including all interaction processes (QEL, RES, DIS, NC). 
In additions, neutrino interactions have been generated in iron to explore the capa- 
bility of the Spectrometer to reconstruct self-contained events taking advantage of its 
additional mass. 

The propagation in the detector is implemented with either GEANT3 or FLUKA 
and the geometry of the detectors is described using the ROOT geometry package. 
The main features of the geometry implemented in the simulation are briefly described 
below. 

The ND LAr target has dimensions 352.8 x 362.5 x 910 cm 3 and a total mass of 
162 t. The target is surrounded by about 1 m thick light material standing for the 
dead material of the cryostat. The basic Spectrometer is an instrumented dipolar 
magnet made of two magnetized iron walls producing a field of 1.5 T intensity in 
the tracking region; field lines are vertical and have opposite directions in the two 
walls whereas track bending occurs on the horizontal plane. The thickness of the iron 
planes is at present envisaged to be 5 cm. Planes of bakelite RPC's are interleaved 
with the iron slabs of each wall to measure the range of stopping particles and to 
track penetrating muons. The Spectrometer is equipped with 20 planes of 3 rows, 
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each consisting of 2 RPC's. At present, additional high precision Drift Tube detectors 
are not simulated. 

The FD target consists of 2 LAr volumes each with dimensions 352.8 x 362.5 x 
1820 cm 3 . The total LAr mass is 648 t. The FD Spectrometer is assumed to be 
similar to the one for the ND (20 planes each with 5 rows and 3 columns of RPC's). 

The response of the Liquid Argon detector has been sketched by sampling the 
position of the charged particle with 0.5 cm resolution. For the RPC's we assume 
digital read-out using 2.5 cm strip width and a position resolution of about 1 cm. 

6.2 Reconstruction 

A framework based on standard tools (ROOT, C++) has been developed for the 
reconstruction in both the Near and Far detectors. 

The reference frame is defined to have the Z-axis along the beam direction, Y 
perpendicular to the floor pointing upwards and X completing a right-handed frame. 
Event reconstruction is performed separately in the two projected views XZ (bending 
plane) and YZ (vertical plane). 

A simple track model is adopted to describe the shape of the trajectory of tracks 
travelling through the detector. The model is based on the standard choice of param- 
eters used in forward geometry (i.e. intercepts, slopes, particle momentum, particle 
charge, ... ). The reconstruction strategy is optimized to follow a single long track 
(the muon escaping from the neutrino-interaction region) along the Z-axis. 

The reconstruction is performed in the usual two steps: Pattern Recognition 
(Track finding) and Track Fitting. 

The task of the Pattern Recognition is to group hits into tracks. Taking into 
account that most of neutrino interactions generated in the LAr target have just 
one track reaching the Spectrometer (the muon) and assuming a read-out capable of 
avoiding event overlap, we postponed the development of a dedicated algorithm For 
the time being all the hits of the events are associated to a single track. 

The Track Fitting has to compute the best possible estimate of the track param- 
eters according to the track model. A parabolic fit is performed in the XZ plane 
(bending) whereas a linear fit is used for YZ plane (vertical). Particle charge and 
momentum are determined from the track sagitta measured in the bending plane; the 
track fit is corrected for material interactions (Multiple Scattering and energy loss). 
Each spectrometer arm provides an independent measurement of charge /momentum. 
A better estimation of momentum is obtained by range for muons stopping inside the 
Spectrometer. The implementation of a track fitting algorithm based on a Kalman 
filter is eventually foreseen. 



43 



7 Physics performances 



In this section we address the physics performance of the FD and ND Spectrometers 
with and without the B field in air option. 

We studied the performance of the Spectrometer using the Monte Carlo simulation 
detailed in the previous section. In particular we report on the negative-focussing 
option which is the most promising for our purposes. We have generated 10 6 v ^ 
(NC+CC) events and 10 6 (NC+CC) events according to spectra shown in the 



right panel of Fig. 14 Neutrino vertices were generated within the LAr volume and 
each outgoing particle from the vertex was propagated in the LAr + Spectrometer 
geometry. Each event was classified into one of the following category: 

1. Fully contained events. These are events in which the neutrino vertex is 
contained the LAr fiducial volume and the muon track stops in the LAr fiducial 
volume. We considered events with a muon track escaping from the vertex at 
least 20 cm long. 

2. Partially contained in LAr not reaching the Spectrometer. These are 
events with the neutrino vertex within the LAr fiducial volume and the muon 
track which escapes from the LAr fiducial volume which do not intercept the 
Spectrometer. In this case we require a muon track at least 200 cm long in 
order to reconstruct the muon momentum with a good momentum resolution 
by MCS in LAr. 

3. Partially contained in LAr reaching the Spectrometer. These are events 
with the neutrino vertex within the LAr fiducial volume and the muon track 
which escapes from the LAr fiducial volume which do intercept the Spectrom- 
eter. In this case we accept also muon track shorter than 200 cm in LAr since 
the muon momentum information are recovered with the Spectrometer. 

The last sample is further subdivided into events with a muons stopping in the 
Spectrometer and events passing through the whole Spectrometer (or escaping from 
the side). In the first case the muon momentum is reconstructed by range, in the 
second case by magnetic bending. The muon charge sign is always reconstructed by 
magnetic bending, in the air field or in the magnetized iron according to the muon 
energy. 

Figg. 28 and [29]report the fraction of each event topology in the negative-focussing 



option in the Near and Far detector, respectively. 



In Fig. 30 we show the Spectrometer performance in terms of momentum recon- 
struction. In the figure the separate contributions of stopping muon, passing through 
going muons and low energy muons are displayed. The total muon momentum at the 
neutrino vertex is computed adding the contribution of the energy lost in LAr. 



In Fig. [31] we show the Spectrometer performance in terms of charge sign recon- 
struction. In the figure the separate contribution of muon bending in the air magnet 
and in the magnetized iron are shown. 
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Figure 28: Fraction of events collected at the Near site by the LAr and the Spectrometer detetctors. 
The different topologies are detailed in the text. 
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Figure 29: Fraction of events collected at the Near site by the LAr and the Spectrometer deetctors. 
The different topologies are detailed in the text. 
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Figure 30: The obtained momentum resolution percentage including all the selection, efficiency and 
reconstruction procedures as described in the text. The blue dots correspond to the measure performed 
by the magnetic field in air, the red (black) dots correspond to what can be achieved by using the 
slope algorithm in the magnetic field in iron with the two (one) arms. The achievable resolution 
with the dE/dX range corresponds to the green dots. 




Figure 31: The obtained charge mis-identification percentage including all the selection, efficiency 
and reconstruction procedures as described in the text. The blue dots correspond to the measure 
performed by the magnetic field in air, the red (black) dots correspond to the magnetic field in iron 
with the two (one) arms. 
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Once the muon momentum is computed, the neutrino energy is estimated accord- 
ing to the event topology. For QE events the two-body kinematics allow a precise 
reconstruction of the incident neutrino energy from the measured momentum and 
direction of the outgoing muon: 

wec = 1 {Ml - ml) + 2E fl (M n -V)- (M n - 

where M p , M n , m M , E^, p M , cos 6^ V are the proton, neutron, muon masses, muon 
energy, momentum and angle with respect to the incoming neutrino direction and the 
LAr nuclear potential, set at 27 MeV, respectively. For non-QE events the neutrino 
energy is reconstructed adding up the muon energy to the energy of hadronic system, 
with a gaussian smearing a E lE = 0.3 /\/E 

A two-flavor sensitivity plot was computed assuming disappearance and a A% 2 
approach. We compared the FD measured neutrino spectrum with the non-oscillated 
spectrum derived from the ND data. We assumed a 5% overall systematic uncertain- 



ties on the absolute normalization derived from the ND spectrum. Fig. [32] shows the 
90% CL. sensitivity curve for the negative-focussing option assuming 3.5xl0 6 pot 
with the iron magnet (the achievable exclusion limits by including also the magnet 
in air, are in progress). 

The study of the possible physics performances of the Spectrometers with respect 
to the neutrino interactions inside the iron is in progress. 



8 Mechanical Structure 

The two OPERA iron dipole magnets [45J can be taken as an example for the design 
of NESSiE Spectrometers. They are made of two vertical walls with rectangular cross 



section and of the top and bottom flux return yokes, as shown in Fig. 33 Each wall 
is composed of 12 layers of iron slabs 5 cm thick, separated by 11 gaps, 2 cm thick, 
hosting RPC detectors. 

The iron layers are made of 7 vertical slabs, for a total area of (8.75 x 8.2)m 2 . 
Including the top and bottom return yokes, the total height of the magnet is about 
10 m and its length along the beam 2.85 m. The magnet weight is around 1 Kton. 

The slabs, the top and the bottom yokes are held together by means of screws, 
while more screws (about 1/m 2 ) are used to keep slabs straight with spacers to ensure 
the thickness uniformity of the gaps hosting the detectors. 

The Spectrometers are magnetized by coils located at the top and bottom return 



yokes, as shown in Fig. 33 



The installation of each magnet was performed according to the following time 
sequence: 

• bottom yoke installation; 
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Figure 32: The sensitivity plot (t 90% C.L.) for the negative-focussing option assuming 3.5x10® 
pot with and without the air magnet option. Black (red) line: V ^ (v ^) exclusion limit. Note that the 
contribution of the magnetic field in air to the sensitivity has not being included, yet. Blue (green) 
line: old (recent) exclusion limits on from previous (CDHS) and recent (!4-3l)measurements.The 
two filled areas correspond to the present exclusion limits on the from CCFR ( I44!i an d Mini- 
BooNE \W$ experiments (at 90% C.L.). 
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• internal support structure construction; 

• iron/RPC layers installation (one plane in each arm at the time in order to keep 
the structure balanced); 

• top yoke installation; 

• removal of the internal structure. 




Figure 33: OPERA magnet scheme. 



For the NESSIE experiment a very similar setup should be used to profit of the 
available detectors which have been designed to obtain the maximum acceptance 
coverage and the strip signal configuration. The different geometrical requirements 
for NESSiE corresponds to an equal width (7 vertical slabs) and a smaller total 
height (5 rows of RPC's instead of 7) for the Far Spectrometer. The ND will instead 
be constructed with 5 vertical slabs of 3.4 m high. For the OPERA experiment a 
total of 336 slabs, (1.25*8.2) m 2 wide, were employed. Such slabs can be cut into the 



49 



shape needed for the present baseline NESSiE configuration, but a new production 
will also be needed. The top and bottom returns must be redesigned and the copper 
coils modified accordingly. 

9 Magnet Power Supplies and Slow Controls 

The following description is based on the past experiences with the OPERA iron 
magnets (see and references therein). Similar applications are foreseen for the 
present Proposal, with minor changes due to the magnetic field in air not present in 
OPERA and the thicker magnet. 

9.1 Power Supply features 

The magnetomotive force to produce the B field is provided by DC power supplies, lo- 
cated on the top of the magnet. They are single-quadrant AC— >DC devices providing 
a maximum current of 1700 A and a maximum voltage of 20 V. As a single-quadrant 
power supply can not change continuously the sign of the voltage, the sign of the 
current is reversed by ramping down the power supply and inverting the load polarity 
through a motorized breaker. The power supplies are connected to the driving coil 
wound in the return yokes of the magnet by means of short flexible cables. 

Ancillary systems 

• Coil: The coil is made of 100 x 20mm 2 copper (type Cu-ETP UNI 5649-711) 
bars. The segments are connected through bolts after polishing and gold-plating 
of the contact surface. Each coil has 20 turns in the upper return yoke connected 
in series to 20 more turns in the bottom yoke. The two halves are linked by 
vertical bars running along the arm. Rexilon supports provide spacing and 
insulation of the bars. 

• Water cooling: Water heat exchangers are positioned between these supports 
and the bars while the vertical sections of the coil are surrounded by protective 
plates to avoid accidental contacts. More than 160 junctions have been made 
for each coil and the quality of such contacts was tested measuring the overall 
coil resistance during mounting. Cooling ensures an operating temperature of 
the RPC detectors lower than 20° C. 

Current status of the two OPERA power supplies 

Even if the l rst power supply stops with a monthly frequencjj^ the 2 nd power 
supply is instead working without troubles. Altogether the overall downtime of the 
two OPERA magnets is about 0.1%. 

'No firm conclusions about the cause of the failures has been reached after several years of 
investigations. 
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A similar framework is foreseen for the NESSiE proposal. Both the Far and the 
Near magnets will need a single power supply each. 

9.2 Monitored quantities for every magnet 

Here is the list of the foreseen values to be monitored: 

• continuous measurement of the magnetic field strength by Hall probes or pickup 
coils at magnet ramp down (not routinely done in OPERA) 

• electrical quantities: 

1. current: check for maximum / minimum range 

2. voltage: check for maximum / minimum range 

3. ground leakage current: check for maximum range (not automatically done 
in OPERA, via web cam or onsite check) 

• temperatures and cooling: 

1. coil temperatures: check for maximum range 

2. cooling water input temperature: check for maximum / minimum range 

3. cooling water output temperature: check for maximum / minimum range 

4. cooling water pressure, cooling water flux: check for maximum / minimum 
values (not automatically done in OPERA, via web cam or onsite check). 

9.3 Slow Control 

The slow control system of the Spectrometer will be developed to master all the 
hardware related: the magnet power supplies, the active detectors and all the ancillary 
systems. 

According to the experience acquired for the OPERA experiment jl6] the system is 
organized in tasks and data structures developed to acquire in short time the status of 
the detector parameters which are important for a safe and optimal detector running. 

The slow control should provide a set of tools which automatize specific detector 
operation (for instance ramping up of the detector High Voltage before the start of 
a physics run) and lets people on shift control the different components during data 
taking. 

Finally, the slow control has to generate alarm messages in case of a component 
failure and react promptly, without human intervention, to preserve the detector from 
possible damages. As an example the RPC High Voltages have to be ramped down 
in case of any failure of the gas system. 

A possible structure of the slow control can be organized as follows: 
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• one databases is the heart of the system; it is used to store both the slow control 
data and the detector configurations. 

• the acquisition task is performed by a pool of clients, each serving a dedicated 
hardware component. The clients are distributed on various Linux machines 
and store the acquired data on the database. 

• the hardware settings are stored in the database and served through a dynamic 
web server to all the clients as XML files. A configuration manager gives the 
possibility to view and modify the hardware settings through a Web interface. 

• a supervisor process, the Alarm Manager, which retrieves fresh data from the 
database, and is able to generate warnings or error messages in case of detector 
malfunctioning. 

• the system is integrated by a Web Server which allows controlling the global 
status of data taking, the status of the various components, and to view the 
latest alarms. 



10 Detectors for the Iron Magnets 

The Near and Far Spectrometers of NESSiE will be instrumented with large area 
detectors, ND and FD respectively, for precision tracking of muon paths allowing for 



high momentum resolution and charge identification capability. According to Fig. 26 
a resolution of the order of 1 mm is needed for the low momentum muons crossing 
the magnetic field in air. A resolution of about 1 cm is instead enough for the large 
mass detectors that we plan to use within the iron slabs of the Spectrometers. In this 
Section we focus on the latter whilst several different options will be described in the 
next Section for the high resolution detectors which have to be used in the "air" part. 

Suitable active detectors for the ND and FD Iron Spectrometers are the RPC - 
gas detectors widely used in high energy and astroparticle experiments [UJ - because: 

• they can cover large areas; 

• are relatively simple detectors in terms of construction, flexibility in operation 
and use; 

• their cost is cheaper than other other large area tracking systems; 

• they have excellent time resolution; 

• and large counting rate power (in specific operational modes). 

Furthermore a large part of the units to be used are already available. In fact by 
considering the remainders of the OPERA RPC production, about 1500 m 2 of RPC's 
can be used. 
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10.1 RPC's detectors 



We plan to use standard bakelite RPC : two electrodes made of 2 mm plastic laminate 
kept 2 mm apart by polycarbonate spacers, of 1 cm diameter, in a 10 cm lattice 
configuration. The electrodes will have high volume resistivity (10 11 — 10 13 f2 cm). 
Double coordinate read-out is obtained by copper strip panels. The strip pitch can 
be between 2 and 3.5 cm in order to limit the number of read-out channels over a 
large area. An optimization of the strip size and orientation (horizontal, vertical and 
tilted ones) is required in order to define the best performances of the Spectrometer 
(track reconstruction resolution and reduction of ghost hits). RPC's are commonly 
used in streamer mode operation with a digital read-out as described in the following. 

10.2 Ancillary systems 

For the operation of the RPC's in the Spectrometer, ancillary systems are needed, 
namely: 

• a Gas distribution system; 

• a High Voltage system with current monitoring performed by dedicated nano- 
amperometers designed by the Electronics Workshop of LNF. 

• monitoring of several environmental/operational parameters (RPC tempera- 
tures, gas pressure and relative humidity). 



10.3 The Gas system 

Since the overall rate (either correlated or uncorrelated) is estimated to be low (see 



Sect. 12), standard gas mixtures for streamer operation can be used, like for instance 
the one employed in the OPERA RPC's, composed of Ar/tetrafTuorethane/isobutane 
and sulfur-hexafluouride in the volume ratios: 75.4/20/4/0.6. 

Different gas mixtures can be further investigated and an optimization is advisable 
in view of the DAQ system adopted (digital versus analog read-out) and the safety 
regulations in the experimental halls. OPERA RPC's are flushed with an open flow 
system at 1500 1/h corresponding to 5 refills/day. The installation of a recirculating 
system could also be considered, if the gas flow has to be increased to prevent detector 
aeine. 



10.4 The Tracking Detectors for the Near and Far Spectrom- 
eters 

The design and evaluation of the requirements for the two detectors are base on the 
RPC's developed for the OPERA experiment at LNGS. Each of its unit of RPC owns 
dimensions of 2904 x 1128 mm 2 . 
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The Near Spectrometer will have a magnetic field in air followed by a calorimeter 
of interleaved planes of iron and RPC In the ND RPC's will be arranged in planes 
of 2 columns and 3 rows for a total exposed surface of about 20 m 2 . A total of 40 
planes will instrument the Spectrometer and a total of 720 m 2 of detectors are thus 
required. Taking into account the RPC dimensions, this number corresponds to 240 
units. 

The Far Detector design will be very similar to the one conceived one for the 
magnetized Iron Spectrometers of the OPERA experiment. The detector will consist 
of units arranged in 3 columns and 5 rows to form planes of about 50 m 2 . As for the 
Far magnet 40 planes of detectors will be interleaved with iron absorbers for a total 
of about 1800 m 2 (600 units). 

10.5 Production and QC tests 

OPERA RPC's, before the installation, underwent a full chain of quality control tests 
serialized in the following steps: 

• mechanical tests (gap gas tightness and spacers adhesion); 

• electrical tests; 

• efficiency measurement with cosmic rays and intrinsic noise determination. 

The setup, still partially available at Gran Sasso INFN Laboratories, was able to 
validate about 100 m 2 of RPC's per week. 

10.6 Costs 

Plastic laminate can be produced in Italy by Pulicelli, a company located near Pavia 
(the company is currently producing material for the CMS RPC upgrade system). 
The cost of plastic laminate is about 30 %jm 2 . The RPC chamber assembly can be 
done in Italy by the renovated General Tecnica company with an estimated cost of 
about 300 €/m 2 . The overall cost of a complete new RPC production is foreseen in 
about 1.5 M€. 

As for the read-out strips, the costs for 5000 m 2 is around 500 K€. 



11 Detectors for the Air Magnet 



The air-core magnet (Sec. 5.2 ) will be used for the charge identification of low momen- 
tum muons which requires precise measurements of the muon path. A (4+4)-layers 



tracker was simulated (see left plot in Fig. 34) assuming a 0.25 T magnetic field of in 



between. Track reconstruction is required only in the bending plane (XZ). The iden- 
tification of the muon charge is optimized looking at the change of the track slope 
before and after the magnetic field. Different detector resolutions were simulated (see 
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right plot in Fig. 34) and a resolution of the order of 1 mm results as the proper 
choice. 

Several detector options for such precise measurements are briefly discussed in the 
following. More detailed MC simulations and eventual test-beams will be crucial to 
test the capability of the different detectors to reach the required performances. The 
final choice on these precision trackers will depend also on the cost, on the reliability 
and on the possibility of re-using parts from previous experiments. 




Figure 34: Left: Display of a simulated event with momentum p — 500 MeV/c in the air-core 
magnet. Right: the charge mis-identification percentage as a function of the muon momentum 
for different detector resolutions (incoming muons perpendicular to the detector plane). The slope 
difference was used to measure the charge. 



11.1 Drift tubes (OPERA-like) 

In the OPERA experiment the High Precision Trackers (HPT's) are mainly devoted to 
muon identification and charge measurement. HPT's are drift chambers (aluminum 
tubes each with a central sense wire) arranged in fourfold layers in order to avoid 
ambiguities in the track reconstruction and to enhance the acceptance (see example 



in Fig. 35 ). In order to simplify the calibration procedures the wires have been located 
by using cover plates and the wire position is decoupled from the position of the tube. 
Thus a wire position accuracy has been achieved with a tolerance of 0.1 mm. In this 
configuration and assuming negligible inefficiency, the spatial single tube resolution 
(rms) has been measured to be better than 0.3 mm. 

The main features of these gas devices are mechanical robustness, absence of glue 
to retain the gas quality, signal quality guaranteed by a Faraday cage. In the OPERA 
detector the HPT's are grouped in modules and the layers are staggered to optimize 
the acceptance and to minimize the left/right ambiguities. The tubes are filled with 
a gas mixture (80% Ar, 20% C0 2 ) and run at a pressure of 1005 mbar. TDC units 
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with a Least Significant Bit of 1.5 ns are used to measure the drift time (the time 
spectrum ranges up to 1.6 fis and TDC's cover a double range). 

The large size (8x8 m 2 ) has been the main challenge in designing the OPERA 
detector, with vertical tubes of ~ 8 m. The HPT size is not a problem in designing 
the NESSiE vertical precise trackers because 3 m and 5 m (vertical dimensions) are 
enough for Near and Far detectors, respectively, and the reduced size could allow the 
complete re-use of the OPERA HPT's. 



Figure 35: Four layers of staggered drift tubes allow reconstructing the track without ambiguities. 



11.2 RPC's with analog read-out 

Resistive Plate Chambers (RPC's) are gas detectors [47J widely used in high energy 
and astroparticle experiments. A single gas-filled gap delimited by bakelite resistive 
electrodes is the simplest set-up commonly used in streamer mode and digital read- 



out. We already introduced them in Sect. 10.1 The most relevant features are the 



excellent time resolution and the high rate capability. Also the position resolution 
is very good, in particular conditions the centroid of the induced charge profile was 
determined [48J with a FWHM resolution of ~ 0.12 mm. In the case of the NESSiE 
experiment such resolution is not necessary and simpler and cheaper set-up can be 
used. 

The analog read-out of RPC has been implemented in the last years [19] or vari- 
ously proposec^] With this technique by reading the total amount of charge induced 
on the strips detailed information can be obtained on the streamer charge distribution 
across the strips and better estimate of the track across the detector is thus achieved 
than in the digital case. The charge profile can be approximated by a Gaussian shape 
whose width (~ 5 mm) does not depend on the gas mixture and operating high 



voltage, unlike the total charge which is strongly dependent on them (see Fig. 36). 

JThe analog read-out of RPC's strips was for example proposed as an alternative option for the 
Target Trackers of the OPERA experiment [50] , 



56 



4000 - 



o 



c 

> 
CD 

o 





1500 h 





49% Ar + 7% iC 4 H 10 


+ 40% C2H 2 F 4 + 4% SFg 




- o 


70% Ar + 20% 


C 4 H 10 + 10% C2H 2 F 4 




fa) 














- 




.-■°* 








- 


p'' 












P'" 

'-X-* — r- -*"T 


— *~ 


• • - 






• 


HV = 8600 V 




I 




* 1 ^ 

(b) 


* 


HV= 9000 V 




„ © 




o 


HV = 9600 V 


,-°* 


'b... 








y 




. 








*■''" 






1 


i'jo. 




10 -5 







5 


10 



mm 



Figure 36: Charge profile (normalized per event) observed with RPC working in streamer mode: 
(a) for various gas mixtures at nominal running HV, (b) for different high voltages with the gas 
mixture (49% Ar + 7%iC 4 H 10 + 40%C 2 H 2 F 4 + 4%SF 6 ) foreseen for the ALICE detector. The curves 
are Gaussian fits of the data. The standard deviation of these distributions is about 5 mm in all 
cases U)l\l. 



A few mm resolution in the charge position determination is obtained by choosing 
an adequate strip size. Also the dynamic range is improved, allowing the detection 
of particles at a density of the order of 1000 particles / m 2 . 

11.3 RPC's with digital read-out in avalanche regime 

Another possibility is to operate the RPC's in avalanche mode: the incoming particles 
release primary charges followed by Townsend avalanches in the gas gap. 

The bidimensional measurement of the avalanche position with mm-accuracy has 
been verified in single-gap chambers [52] . As in avalanche regime the charge cluster 
size results to be smaller by some order of magnitude with respect to the streamer 
mode operation (it depends on HV and gas mixture), it is possible to achieve better 
space resolution by using smaller digital strips. A higher rate capability is also at- 
tainable due to the lower amount of charge delivered in the avalanche. The response 
of standard 2 mm gap RPC's in streamer mode can not be higher than 3 MHz/m 2 , 
while in avalanche mode it can attain 30 MHz/m 2 . 

It is also possible to use the RPC in saturated avalanche mode. The advantage is 
to have a lower charge signal than in streamer mode but still high enough to remove 
the need for a pre-amplifier. Analog read-out with an optimized strip size guarantees 
an adequate position resolution. 
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Anyway the avalanche mode has the inconvenient of being more sensitive to tem- 
perature and pressure variability due to the lower amount of charge produced. 

Finally, regarding the possibility to get high precision position measurement with 
RPC's at a sustainable cost, a new procedure has been developed to determine the 
charge position by timing measurements. Precision in the sub-millimeter range has 
been reached [53]. The method is based on the read out of the signal propagating in 
the graphite. The graphite electrode, coupled to the ground reference of the detector 
read-out panels, is a distributed capacitance-resistance system (see Fig. 37). The 



occurrence of a discharge in the gas produces a point-like perturbation of the steady 
potential distribution of the system. The time behaviour of the perturbation can 
be described, in the approximation of infinite electrode size, by a two-dimensional 
Gaussian distribution with time dependent variance and amplitude. The distance is 
related to the time (of the maximum) by a quadratic relation which can be used to 
measure the position from the time measurement. The sharpness of the maximum 
is strongly dependent on the distance. A distribution with a FWHM of 0.8 mm has 
been reached. 




Figure 37: Pick-up scheme in the graphite electrode. From 



12 Backgrounds 

Assuming neutrino fluxes as described in Sect. [4] around 10-15 events per spill are 
expected in the Near Detector, a spill being 2.1 /xs long within a cycle of 1.2 sec. For 
the Far Detector the rate of events is reduced by about a factor of 20. 

The possible background rates are analyzed assuming a tracking system with 
single-gap RPC's. We distinguish the uncorrelated background due to detector noise 
and local radioactivity (dark counting rate) and the correlated background due to 
cosmic rays. 
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12.1 Uncorrelated background 

The dark counting rate depends on detector features and ambient radioactivity. A 
typical value for RPC as measured at sea level is 300 Hz/m 2 . Therefore the expected 
rate per plane is \n = 300 Hz/m 2 x (3 x 6 m 2 ) ~ 5.4 kHz on the Near Spectrometer 
and Xp = 300 Hz/m 2 x (5 x 9 m 2 ) ~ 13.5 kHz on the Far Spectrometer. Assuming 
a read-out time window of 10 fis we expect 0.05 (0.14) hits per RPC plane per event 
on the Near (Far) detector. The number of fired strips will depend on the strip width 
(in OPERA with 2.6 and 3.5 cm strip-wide the typical cluster size is ~ 1.5 strips). 

The requirement of three-contiguous-planes coincidence in the beam-spill time 
(tipically 2 fis) makes the dark noise contribution to the trigger rate negligible (see 
last column in Tab. [8]). 

12.2 Cosmic Ray background 

The contribution of Cosmic Rays (CR) to the plane-by-plane background is similar to 
that discussed in the previous Section, but CR events yield long tracks that constitute 
a potentially more dangerous correlated background. 

Assuming a trigger majority of at least 3 fired planes a cosmic particle can trigger 
the data-taking when it has enough energy to cross 2 iron slabs (2x5 cm) . The CR flux 
at sea level is essentially due to muons (hard component of Extensive Air Showers) and 
these particles can trigger the data-taking when their momentum exceeds 250 MeV/ c. 
Then the integrated vertical flux of muons is J = 97 Hz/m 2 / sr at sea level [54J. 

The integrated vertical flux of the soft component (electrons and positrons) can 
be represented by J(> E) = 0.22 E~ 1A5 [Hz/m 2 / sr], where E is in GeV. Taking 
into account this contribution and minor ones due to hadrons, J = 100 Hz/m 2 j sr 
will be used in the following calculations as integrated vertical CR flux at sea level. 

In the conservative hypothesis that the CR ray flux is isotropic above the horizon 
and equal to the vertical flux, the total rate \rc on a detector shaped as a fully 
efficient box is 

Xrc = \S W J 

where S to t is the surface of the detector. The expected number of CR events in a time 
window of 2 fxs (the beam-spill time) is reported in Tab. [8] They scale with the time 
window T (by a factor T/2 fis). The data in Tab. [8] conservatively ignore that more 
detailed trigger conditions allow a significant reduction of the background. 
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Dimensions Surface Total Iron Mass 
(m x m x m) (m 2 ) (ton) 


CR events Dark noise ev. 
in 2 fis in 2 fis 


Near 
Far 


6 x 3 x 5 126 364 
9 x 5 x 5 230 846 


4 x 10"2 2 x 10" 4 
7 x 10" 2 3 x 10" 3 



Table 8: Dimensions and surfaces of the detectors, iron masses and estimated background events 
(CR and dark noise) for the Near and Far Spectrometers. 



13 Read-out, Trigger and DAQ 
13.1 DAQ overview 

The aim of the DAQ is to read the signals produced by the electronic detectors and 
to create a database of detected events. We recall in Table [9] the characteristics of 
the CERN PS primary proton beam. These are important inputs to define properly 
the data acquisition and flow. 

Table 9: Characteristics of the PS primary proton beam. 





PS Parasitic 


PS Dedicated 


Proton beam momentum 


20 GeV 


20 GeV 


Protons per pulse 


2.6 x 10 13 


3 x 10 13 


Number of bunches 


7 


8 


Bunch length (4 sigmas) 


65 ns 


65 ns 


Bunch spacing 


262 ns 


262 ns 


Burst length 


1.8 fis 


2.1 fis 


Maximum repetition rate 


1.2 s 


1.2 s 


Beam energy 


84 kJ 


96 kJ 


Average beam power 


70 kW 


80 kW 



The foreseen DAQ architecture is composed of three stages: 

• the front end electronics close to the detector (FEB) 

• the read-out interface which together with the trigger board control the read-out 
of the FEB 

• the Event Building which reconstructs events using standard workstations. 

The whole event reconstruction is based on the time correlation of the channels, which 
depends on the accuracy of the electronic channels time stamping. A time resolution 
in the range 5 -j- 10 ns is sufficient to correctly associate the different hits to the 
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corresponding events. A common time reference with respect to the PS extraction 
time is used to time-stamp the data and to correlate the events recorded by the 
Spectrometer to the data recorded in LAr. 



13.2 Data Flow 

The Far detector is designed with 15 RPC/plane x 20 planes/arm times 2 arms for 
a total of 600 RPC detectors, 3 m 2 each. The Near detector is designed with 6 
RPC/plane x 20 planes/arm times 2 arms for a total of 240 RPC's. 

Assuming a strip size of 2.6 cm along Y and 3.5 cm along X for the Far detector 
each plane will be equipped with 32 x 5 = 160 horizontal strips plus 112 x 3 = 336 
vertical strips for a total number of about 500 electronic channels per plane. The total 
number of channels is therefore about 20,000. For the Near detector the number of 
electronic channels will be about 300 per plane. The total number of channels is 
about 12,000. 

Given the expected background rate due to the RPC single rate and the counts 



due to cosmic ray within the beam spill reported in Sect. [L2j the expected data rate 
is dominated by the beam related events. With a maximum PS beam intensity of 
3 x 10 13 p.o.t. about 10 to 20 events are expected in the Near detector. Assuming 
four hit strips per plane and at most 16 byte of data per hit (channel address, signal 
and time) the size of the event after zero suppression is expected to be 1 kbyte. 



13.3 Front-End Electronics 

The role of the electronic read-out is to discriminate the signals coming from the 
RPC's strips and to record the signals above threshold. The start-of-burst signal will 
be used as a trigger and discriminated signals will thus be recorded during the whole 
burst duration. A time-stamp with a resolution of 5 -j- 10 ns will allow associating 
hits belonging to the same events. 

The electronic read-out of the RPC's could be developed according to the same de- 
sign scheme adopted for the OPERA experiment, where they are operated in streamer 
mode and read-out by means of twisted pair cables; these are in turn connected to 
special interconnection boards allowing the front-end electronics to be placed far from 
the detector. The Front-end Boards (FEB) collect signals from RPC's strips and de- 
liver them to a Controller Board (CB) acting as an interface to the DAQ system. 

About 32,000 channels will be used to read-out 2,500 m 2 of RPC detectors with 
pickup strips of 2.5 cm pitch in the vertical direction (orthogonal to the bending plane) 
and 3.5 cm pitch in the horizontal direction (tracking without bending)]^] Each FEB 
board will collect signals from 64 strips providing discrimination, time stamping and 

k To get a spatial resolution of about 2 mm strips with a pitch of 1 cm have to be used for the 
measurement in the bending plane. In such case the total number of channels increases of about a 
factor 2. 
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local buffering sized to store up to 40 events per burst. The local buffer will be read- 
out from the CB, during the inter-spill time, lasting 1.2 s. Serial FEB data transfer 
to the CB will occur at a clock of 10 MHz, taking a few hundred microseconds. A 
new ASIC front-end chip will be designed to assign the time-stamp to each event and 
make it possible to record multiple events in the same burst. The possibility to have 
a double electronic chain, digital and analog, for each channel will be considered. 
In the ASIC project a 10-bit ADC with a multiplexing and sample-and-hold system 
allowing for the analog read-out could be added. The design of this new ASIC chip 
can be provided by the INFN electronics CAD service. 

The cost estimation, based on a production of 2,500 16-channel chips, amounts to 
30 €!/chip. The production cost of 700 Front-End Boards, including additional com- 
ponents and assembly, has been estimated to be about 50 €/board. The production 
of 100 VME Controller-Board has been estimated to be 300 €/board. Additional 
costs: cables and connectors 20 K€, low voltage power supply 20 K€, crates 30 K€. 
The total cost of the RPC's read-out electronics has been estimated to be 210 K€, 
taking into account 32000 digital channels plus some contingency. 

13.4 DAQ 

The acquisition for each detector (FD/ND) is composed by Front-End Boards (about 
5 FEB per plane), Trigger Boards, Controller Boards (1 CB per plane) and the Event 
Builder. In order to acquire with a rate of 10 to 20 events/spill (each spill is 2.1 
fis long), Front-End Boards should provide a pipeline system to save all the data of 
a spill, provide a FAST.OR signal for trigger purposes of granularity in the range 
between 8 and 32 strips and perform the zero suppression. The possibility of a digital 
or analog read-out must be envisaged. The start of the pipeline should be given by 
the start of the spill and the time resolution should be in the range 5-10 ns. 

The Trigger Boards should process the FAST.OR signals from the Front-End 
Boards and provide an external trigger condition to the Front-End Boards. This 
signal could also be used as an external trigger for the High Precision Tracker read- 
out (t signal for the TDC). A Trigger Bus should be implemented for the distribution 
of the trigger. Given the expected overall rate (beam events, cosmic, single rate) and 
due to the Front-End Boards buffer a trigger signal based only on the beam spill 
could be envisaged. 

The Controller Boards should read-out the signals from the Front-End Boards and 
propagate them to the Event Building which is receiving data from the Spectrome- 
ter and the High Precision Tracker systems. The Event Building is using standard 
network protocols. 

Data read-out and Event Building should be performed in the time between two 
spills (~ 1.2s). The Event Builder should be based on standard commercial worksta- 
tion. Data spying and monitoring process will also be implemented at this level. 
The correlation between Spectrometer and LAr data is achieved using a common 
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clock signal to time-stamp the events. Data merging is performed offline at the re- 
construction level. 

14 CERN Logistics 

The two experimental halls identified in the LAr proposal (see Sect. 3 of |lj) to host 
the Near and Far detectors, B181 and B191, respectively, are well suited for the two 
Spectrometers too. The two halls have been chosen to be in the line of the PS neutrino 
beam and to fit the size of the detectors. Some discussion is undergoing with CERN 
for the use of the chosen pit in hall B181, already exploited by some servicing. As 
the pit would not actually fit both the LAr and the Spectrometer detectors and some 
excavation works will be anyhow needed, a different solution is in addition under 
study. A new pit in a new hall in front or behind B181 may be constructed. 




Figure 38: The general layout of the Hall B-191 at CERN that is supposed to host the Far detectors, 
LAr and Spectrometer. The Spectrometer is placed behind the LAr tanks and overlays the cryogenic 
plant, which has to be moved e.g. on the right side. The direction of the neutrino beam is also 
indicated on the top of the picture. 
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For the Far site (see Fig. 38) the Spectrometer detector can easily be allocated 
behind the two Icarus tanks in case the cryogenic equipment be allocated laterally. 
The total weight for the magnet will exceed 1 Kton, on an extension of less than 40 
m 2 . Therefore the floor has to be eventually adapted. As described in Sect. [8] the 
installation of the magnets has to follow a precise time sequence which has to be com- 
patible with the LAr installation and contemporary allow for the closest positioning 
of Spectrometer and LAr. Detailed engineering studies will take care of this subject. 

15 Schedule and Costs 

The choice of performing a design study that is reliable under several aspects led to 
make conservative, well controlled and realistic options. In particular the choice to 
use detectors like the RPC ones and the development of dipole magnets which had, at 
least partly, already been constructed and are in use, allows us to keep under control 
both the time schedule and the costs. 



With respect to the schedule, which is reported in Tab.[10j it should be underlined 
that is based on the deep experience acquired with the OPERA Spectrometers, built 
up from 2005 to 2006 under critical condition^] 



Year(portion) 


Action 


2012 


Design optimization 


l rst half 2013 


Define tenders/contracts 
Setting up Detectors Test-stands 


2 nd half 2013 


Mechanical Structure construction 
Detector production 
Start Detector test 
Magnet (Air) preparation 


Y st half 2014 


Start Magnet installation 
Start detectors installation 


2 nd half 2014 


End installation 



Table 10: Tentative time schedule for the 2012-2014 years which will include the final optimization 
of the detectors, their construction and installation at CERN. 

With respect to the cost estimate, the expenses needed for the major items are 



reported in Tab. 11 



The partial re-use of systems developed for the OPERA experiment may be en- 
visaged in case that experiment would not need them anymore at the time of the 
construction. We note that the OPERA Spectrometers have been fully funded by 



'During the period 2005-2007 LNGS laboratory underwent restoration and safety works enforced 
by Italian Government, following the temporary seal of May 2005. 
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Item Cost (in M €) 



Far 


Magnet 


2.5 


Detectors 


1.0 


Strips 


0.5 


Data Acquisition 


1.0 


Near 


Magnet 


2.0 


Detectors 


0.5 


Strips 


0.3 


Data Acquisition 


0.5 


Total 


8.3 



Table 11: Estimate of the costs of the major items. 

INFN, except the Precision Trackers, which is therefore committed to their disman- 
tling and entitled to possibly re-use as well. To this respect, it is clear that e.g. the 
iron raw material will be available sooner or later to INFN which may then take it 
into account for company arrangements. 

16 Conclusions 

Existing anomalies in the neutrino sector may hint to the existence of one or more 
additional sterile neutrino families. We performed a detailed study of the physics 
case in order to set a Short-Baseline experiment at a refurbished CERN-PS neutrino 
beam, able to either prove or reject the existence of sterile neutrinos. 

The already submitted proposal based on the technology of imaging in ultra- 
pure cryogenic Liquid Argon (LAr) may suffer from some experimental limitations 
which we deem critical: the measurement of the muon charge on event-by-event basis 
extended to the lowest achievable energy range would be mandatory. Indeed the muon 
leptons from Charged Current (CC) (anti)-neutrino interactions play an important 
role in disentangling different phenomenological scenarios, provided their charge state 
is determined. Since CC muon events induced by the primary beam would be very 
abundant, a search based on muon appearance-disappearance will benefit of large 
statistics. 

The best option in terms of physics reach and funding constraints is provided by 
two Spectrometers based on dipoles magnets mostly in iron, at the Near and Far sites 
(located at 127 and 850 m from the CERN PS neutrino beam, respectively), to be 
placed behind the LAr detectors. 

In order to measure the momentum and identify charge with high sensitivity in 
an extended energy range (from few hundred MeV to above 5 GeV) we complement 
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the iron-core dipole by a magnetic field in air in a limited region just in front. 

The selected detectors, mainly Resistive Plate Chambers, would exploit well 
known technologies and make it possible to re-use part of existing ones (should they 
become available; if not, it would imply an increase of the costs with no additional 
delay) . 

The Near and Far Spectrometers by complementing the capabilities of the LAr 
detectors will fully exploit the very rich and exciting physics information in case sterile 
neutrinos indeed exist. In fact the Spectrometers by measuring the momentum and 
identifying the charge of the muons will provide valuable information in: 

• measuring disappearance in the full momentum range which is a key in- 
gredient in rejecting the anomalies or measuring the whole parameter space of 
oscillations involving sterile neutrino; 

• measuring the neutrino flux in the Near detector, in the full muon momentum 
range, which is quite relevant to keep the systematic errors low. 

The measurement of the muon charge will furthermore provide valuable information 
in: 

• separating from in the antineutrino beam where the contamination is 
quite important. This measurement is critical to fully exploit the experimental 
capability of measuring a difference between — > u e and z/ M — > V e difference in 
view of a possible signature of CP violation. 

Results of our study are reported in detail in this proposal where our experiment 
is identified with the acronym NESSiE (Neutrino Experiment with Spectrometers in 
Europe). 
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